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m ' ABSTRACT 

in ; 

i The fluctuations in the Cosmic Microwave Background (CMB) intensity due to the 

' Sunyaev-Zeldovich (SZ) effect are the sum of a thermal and a kinetic contribution. 

Separating the two components to measure the peculiar velocity of galaxy clusters 
• requires radio and microwave observations at three or more frequencies, and knowledge 

Q-f of the temperature T e of the intracluster medium weighted by the electron number 

q density. To quantify the systematics of this procedure, we extract a sample of 117 

;_h ' massive clusters at redshift z = from an iV-body hydrodynamical simulation, with 

\ 2 x 480 3 particles, of a cosmological volume 192 h~ x Mpc on a side of a flat Cold Dark 

Ci • Matter model with ilo = 0.3 and 17a = 0.7. Our simulation includes radiative cooling, 

star formation and the effect of feedback and galactic winds from supernovae. We 
find that (1) our simulated clusters reproduce the observed scaling relations between 
X-ray and SZ properties; (2) bulk flows internal to the intracluster medium affect the 
5— i ' velocity estimate by less than 200 km s _1 in 93 per cent of the cases; (3) using the X- 

ray emission weighted temperature, as an estimate of T e , can overestimate the peculiar 
velocity by 20 — 50 per cent, if the microwave observations do not spatially resolve 
the cluster. For spatially resolved clusters, the assumptions on the spatial distribution 
of the ICM, required to separate the two SZ components, still produce a velocity 
overestimate of 10 — 20 per cent, even with an unbiased measure of T e . Thanks to 
the large size of our cluster samples, these results set a robust lower limit of ~ 200 
km s _1 to the systematic errors that will affect upcoming measures of cluster peculiar 
velocities with the SZ effect. 

Key words: large-scale structure of Universe - galaxies: clusters: general - cosmology: 
miscellaneous - methods: numerical 
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1 INTRODUCTION 

The properties of the spatial distribution and velocity field 
of matter in the Universe constrain the models of large- 
scale structure formation. The velocity field has mostly been 
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used to estimate Qo, the mass density parameter of the 
Universe, by means of the clustering anisotropy in red shift 
space (e.g. iHamiltonlll998l ; iMatsubara fc Szalavl liboa) . the 
connec tion between the densit y and the peculiar velocity 
fields llCourteau fc D ckcl 200j and references therein), or 
the mean relative peculiar v elocity of galaxy pair s as a func- 
tion of pairwise separation jFeldman et al.ll2003l) . 

Further constraints on structure formation models can 
be extracted from the evolution of the velocity field on large 
scales. Its measure is possible with the detection of the ki- 
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netic SZ effect iSunvaev fc Zeldovichlli972fr . When the Cos- 
mic Microwave Background (CMB) photons are scattered 
by the free electrons in the intracluster medium (ICM), the 
CMB radiation varies its intensity and its spectral distri- 
bution. The thermal motion of the free electrons and their 
bulk velocity have a different spectral signature on the CMB 
and the two contributions can be separated with multi- 
frequency observations. However, the bulk velocity contri- 
bution (the kinetic SZ effect) is, on average, an order of 
magnitude smaller than the thermal contribution (the ther- 
mal SZ effect), and leaves the CMB spectrum unchanged: 
these features make the kinetic SZ effect difficult to detect, 
although it has been recently s hown how its spatial correla- 
tion with the thermal SZ e ffect jPiaferio. Sunvaev fc Nusserl 
120001: Ida Silva et al]l200ll) can be used to remov e the CMB 
contribution efficiently jForni fc Aghanimll20o"3) . 

The thermal effect is now becoming routinely detected 
and major efforts have been done to plan several SZ cluster 
surveys which will produce catal ogues of tens of thousands of 
clust ers in the comin g years (see lCarlstr om. Holder fc Reese! 
120021 for a review and lSchulz fc Whitel2003T for systematics) . 
The SZ effect does not suffer from the usual redshift dim- 
ming of electromagnetic surface brightnesses, and clusters 
can be detected, in principle, as far back in time as when 
they formed. 

The distribution of SZ clusters on the sky can 
alread y provide informatio n on cosmological parame- 
ters jMei fc Bartlett 120031. 2004) and cluster biasing 
JPiaferio et alJ 120031: Icohn fc Kadotal 120041) . although one 



has to separate the effects of the ICM properties f rom those 
of the cosmological model jMoscardini et al I2002D . Stronger 
constraints for cosmology are of course provided by the 
clust er red shift distribution (e.g . iHolder^Jaim an fc Mohr l 
| 200ll: iBenson. Reichardt fc Kamionkowskil 120021: 
IWeller. Battve fc KneissJ l2002rl . The accuracy of clus- 
ter redshifts determ ined with photometric observations 
jHuterer et al.ll2004l) or methods based on the SZ informa- 
tion alone foiego^tal]|2003t ISchafer. Pfrommer fc Zaroubil 
2003) may be sufficient to constrain cosmology, and optical 
or X-ray follow-up spectroscopy, which can be particularly 
demanding in terms of observation time and sensitivity, 
may not be needed for a first analysis. 

If we measure the peculiar velocity of clusters, 
we can reconstruct the redshift evolution of the ve- 
locity field on large scales. Many authors have in- 
vestigated the feasibility of measuring the bulk flow 
on ~ lOOfe -1 Mpc scale s with the kinetic SZ effect 
llHaehnelt fc Tegmarklll996l: iKashlinskv fc Atrio-Barandeial 



2000; 
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Atrio- Barandcla. Kashli nskv fc Miicketl 12004) . No other 
method is currently viable for this task other than the 
kinetic SZ effect. 

To date, because of the weakness of the total SZ sig- 
nal compared to the sensitivity of the radio telescopes, sep- 
arating the kinetic from the thermal contribution has re- 
quired the estimate of the expected SZ signal with an a 
priori knowledge of the temperature and spatial distribu- 
tion of the ICM. X-ray observations have been used to ob- 
tain this information for spatially well resolved clusters (e.g. 
[ Holzapfel et alJ Il997l: iLaRoaue et al.l 120031 : IBenson et al] 
120031. 2004: iKitavama et alJl2004h .~ 

This procedure can not be easily extended to clusters 



at high redshift where X-ray observations can be beyond 
the limit of detectability and spatial resolution can be poor. 
In principle, but only for massive clusters with hot ICM 
(> 8 keV), one can estimate the temperature, rather than 
with X-ray observations, by using the dependence of the 
spectral function of the SZ effect on the gas temperature, 
which appears when a fully relativistic treatment of the pho- 
ton sc atter is considered iPoi ntecout eau. Giard fc Barred 
119981 : lltoh. Kohvama fc Nozawal Il998t) . However, a ~ 
cent attempt, although successful, has provided a 
temperature measure with q uite large uncertainties 
jHansen. Pastor fc Semikodl2002T) . 

Further insights into the physical properties of the 
ICM come from the scaling relatio ns between X-ray 
and S Z observables. Analytic m odels llDos Santos fc Dorel 
12001 iMcCarthv et all l2003al). semi-anal ytical models 
JCavaliere fc Mencill200lt IVerde et"aH 120021) . and TV -body 
hydrodynamical simulations fe-g- TdgH^nva et alJl2004l) have 
been used to investigate the effects of heating, cooling, and 
preheating on the normalization, slope and scatter of these 
scaling relations. 

Here, we use a la rge iV-body hydrodynamical simulation 
feorgani et all2 004) , which includes an advanced treatment 
of the gas dynamics, contains a large sample of simulated 
massive clusters and reproduces reasonably well the X-ray 
properties of real clusters, to verify (1) whether we can re- 
produce the observed scaling relations between X-ray and 
SZ observables; (2) whether the two basic assumptions un- 
derlying the peculiar velocity measurement hold in the simu- 
lations; namely whether (i) the mean gas velocity equals the 
dark matt er bulk velocity and the internal gas bulk flows are 
negligible jNagai. Kravtsov fc Kosowskvll2003h . and (ii) the 
ICM temperature derived from X-ray observations is a good 
estimator of the ICM temperature weighted by the electron 
number density. Our simulated sample contains 117 well re- 
solved clusters and provides a major statistical improvement 
on previo us analyses performed on a hand ful number of clus- 
ters fe.g. lNagai et~ai1l2003l:lHoldeill2004l) . 

We do not investigate the systematics due to a real- 
istic observational procedure, which includes cleaning the 
radio detection of instrumental noise, atmospheric emis- 
sion, galactic emission at low frequencies, galactic dust and 
young galaxy emission at high frequencies, radio sources, 
and the CMB itself. A det ailed analysis of these effe c ts can 
be f ound, for example, in [Knox, Holder fc Church! (l200-4h 
and lAghanim. Hansen fc Lagachd <l2004t) . Our results will 
only provide a lower limit to the systematic uncertainties on 
the peculiar velocity measures. 

The layout of this paper is as follows: in Sect.|2]we sum- 
marize the basic equations used to estimate the peculiar ve- 
locity with the SZ effect; in Sect. [2] we describe the A-body 
simulation and our sample of simulated clusters; Sect. ^dis- 
cusses the scaling relation between X-ray and SZ observables 
and Sect.|S|tests the basic assumptions on the velocity mea- 
sures. In Sect. we estimate the systematic errors affecting 
the peculiar velocity measurements. Our conclusions are in 
Sect.|3 
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2 BASICS 

Here, we describe the basic equations governing the SZ ef- 
fects and their observational analysis. The hottest cluster in 
our simulated sample has temperature as 7 keV and we can 
safely ignore any relativistic corrections which are relevant 
for clusters with an ICM temperature > 8 keV. 



2.1 Resolved Sources 

At a given frequency u, the variation in the CMB specific 
intensity due to the thermal motion of free electrons in the 
ICM is, in the limit of non-relativistic electron speeds (e.g. 
lRephaelilll995l) . 



AI t (x,8)=i g(x)y(e) 



(1) 



where x = hpv/kT, k is the Boltzman n constant, T = 2.72 5 
K the present day CMB temperature jMather et al] ll999'). 
hp the Planck constant, i = 2(kT) 3 /(h P cf = 1.129 • 
10 3 (T/Kf mJy arcmin" 2 = 13.34(T/7T) 3 MJy sr _ \ c the 
speed of light, and 



Q\x) — -, TT X 4 



(2) 



the spectral function. The Comptonization parameter is 
y{0) = Mfl)) los fc(reW 2 )los , 
where m e is the electron mass, 



(3) 



(t(0))ios = ctt / n e (r)dl 



(T e (6»))i os 



/ n e {r)T e {r)dl 
J n e (r)dl 



(4) 



(•») 



and (jt is the Thomson cross section; n e and T e are the 
electron number density and temperature, respectively. The 
integrals are over the line of sight I, and r 2 = d 2 A 6 2 + I 2 , 
where <1a is the angular diameter distance to the cluster, and 
is the angular separation vector from the cluster centre. 
(T e (6))\ OB is the optical depth weighted temperature and it 
generally differs from the pressure-weighted temperature 



(T p (6))i OB — 



J n e {r)T 2 {r)dl 
f n e (r)T e (r)dl 



(6) 



{T p (8))\ os can be derived from SZ observations alone by in- 
cludi ng the relativ i stic correction in the spectral function 
g(x) <HanserJl2061 iKnox et al]l2004t) . In the following, we 
will not consider this temperature because our simulated 
clusters are not massive enough to make the relativistic cor- 
rection relevant; therefore, with the sensitivity of the tele- 
scopes available now or in the near future at the radio and 
microwave wavelengths, the pressure-weighted temperature 
(T p (Q))\ os of our simulated clusters will not be measurable. 

The bulk flows of both the ICM within the cluster and 
the cluster itself yield a further variation in the CMB inten- 
sity: 



AI k (x,0) = -i h(x)b(0) 

x 4 e x 

h{x) = W^W ' 

b(0) = <t(0))ios</3(0)) 1os , 



(7) 
(8) 
(9) 



1 f n e (r)v los (r)dl 
c J n e (r)dl 



(10) 



One actually observes the sum of the thermal and ki- 
netic contributions 

AI(x, 6) = i [g(x)y(9) - h(x)b(9)] . (11) 
At 6 = 0, we have 

AI {x) = i [g{x)y - h(x)b ] , (12) 

where y = T (kT e /m e c 2 ), r = (r(0))i os , T e = (T e (0))i os , 
and b = to/3, (3 = (/3(0)) los . 

Usually, to extract the SZ image of the cluster from the 
radio/microwave data, a /3-model of the gas distribution is 
assumed. The shape parameters f3 c and 6 C are derived from 
a fit to the X-ray surface brightness 

2 \ l/2-3/J e 



AS(0) = ASo (1+^2 



(13) 



At a frequency x, the CMB intensity variation due to the 
SZ effects is expected to be 



AI(x,e) = AI {x) 1 + 



e_f_ 

o- 



(14) 



where A/o(a:) is now the only parameter provided by the mi- 
crowave observations. When we have a measure of A7o(a:) at 
three or more frequencies x, a fit to equation l)12[l provides 
an estimate of yo and bo. Simultaneous multi-frequency ob- 
servations can be performed to minimize the variations in 
the atmospheric emission. In this case the fitting procedure 
to the spectral function is obviously considera bly more com- 
plicat ed, because the A7o(s)'s are correlated feenson et alJ 

By combining yo and bo, one finally derives the line-of- 
sight peculiar velocity: 



kTe bo 



(15) 



m e C y 

Provided that the spherical j3— model adequately describes 
the gas distribution, when applied to real clusters, this equa- 
tion has two main problems, that we will investigate below: 
(1) the central value of the velocity j3 = {/3(0))i os might not 
necessarily coincide with the global peculiar velocity of the 
cluster; (2) only the temperature averaged within the region 
of the size of the beam is actually available, rather than the 
temperature in the cluster centre T e = (T e (0))i os . 

2.2 Unresolved sources 

If the spatial resolution of the radio telescope is not sufficient 
to resolve the source, the observable quantity is the variation 
of the CMB specific intensity due to the SZ effect integrated 
over the beam size 9t (we assume here an ideal step function 
window) 

AS(x,e b )= / AI{x)dn = io{g(x)Y(8 b )-h(x)B(e b )} ,(16) 
where 

Y(0 h )= f y(0)dfi=(r) flb ^% 



(17) 



4 Diaferio et al. 



B{O b )= / b(8)AV = (T) Bb (l3) 6b 



(18) 



(T)e b =cTT n e (r)dldn = -f / n e (r)dV, (19) 



and 



X,( 9i ,) n ^ r ) r -( r ) d3r 

X.(« b ) ne(r ) d3r 



f Co N n e(»")d 3 



(20) 



(21) 



where rA(#f,) is the radius of the solid angle circle in physical 
units, dl is the proper coordinate element along the line of 
sight, and dPr — d 2 A dldQ is the proper volume element. With 
radio/microwave observations at three or more frequencies, 
we can determine Y(0b) and B(9 b ) as free fit parameters. 
The peculiar velocity is thus 



(0)o b = 



k(T e ) Bb B{6 h ) 
m e c 2 Y(6 b ) 



(22) 



3 THE MODEL 

3.1 The Simulated Cluster Sample 

iBorgani et alJ l)2004h describe in detail the simulation we 
use here. Briefly, we simulate a cubic volume, 192 h" 1 Mpc 
on a side, of a flat ACDM universe, with matter density 
n = 0.3, Hubble constant H = 100ft km s _1 Mpc -1 , 
h — 0.7, baryon density f2bar = 0.02/i~ 2 and power spec- 
trum normalization ag = 0.8. The density field is sampled 
with 480 3 dark matter particles and an initially equal num- 
ber of gas particles, with masses mnm = 4.6 x 1O 9 /i _1 M 
and m gas = 6.9 x 10 8 /i -1 Mq, respectively. The Plummer- 
equivalent gravitational softening is 7.5 h' 1 kpc comoving 
at z > 2, and fixed in physical units at lower redshift. 

The simulation was r un with GADGET-2 
iSpringel. Yoshida fc White! [2001), a massively parallel 
Tree+SPH code with fully adaptive time-stepping, which 
uses th e energy and entropy conser ving SPH implementa- 
tion of ISpringel fc Hernauist| i2002h . The code includes a 
photoionizing, time-dependent, uniform UV background, 
radiative cooling, star formation, feedback from type II 
supernovae and a phenomenological recipe for galactic 
winds. Within each gas particle of sufficiently high density, 
the gas is a two-phase fluid, with unresolved cold clouds, 
embedded at pressure equilibrium in an ambient hot 
medium, providing a sub-grid mod el for the multiphase 
natur e of the interstellar medium dSpringel fc Hernauisd 
120031) . 

The volume of our simulation is larger than 140 3 /i -3 
Mpc 3 , as required to p redict the evolution of cluster pecu- 
liar velocities correctly JSheth fc Diaferiol200ll) : this volume 
also yields a cluster sample large enough for statistical pur- 
poses. Moreover our simulation represents a substantial im- 
provement compared to recent cosmological hydrodynamical 
simulations used for the study of the SZ effect: they typi- 
cally have a volume W0 3 h~ 3 Mpc 3 and a pa rticle masses 
four to six times larger than in our simulation dWhite et alJ 



Figure 1. Maps along three orthogonal directions of the thermal 
SZ effect for a simulated cluster with mass M(< R v i r ) = 6.95 X 
10 14 h~ 1 M@. Each pixel has size 61.6/i — 1 kpc on a side and the 
field of view is 7.9h~ 1 Mpc on a side. 



Figure 2. Same as Figure HI for the kinetic SZ effect. 



l2002Uda Silva et ahlEooll . Finally, the good mass and spa- 
tial resolution and the treatment of a two-phase fluid for 
star formation and feedback provide a reasonable modelling 
of the gas physics. 

We identify clusters in the simulation box with a two- 
step procedure: a friends-of-friends algorithm applied to the 
dark matter particles alone provides a list of halos whose 
centres are used as input to the spherical ov erdensity algo- 
rithm which outputs the final list of clusters iBorgani et alJ 
2004). Centered on the most bound particle of each cluster, 
the sphere with virial overdensity A e (z), with respect to the 
critical density feke. Cole fc Frenldll996h . defines the virial 
radius R v i T . For this cosmological model A c = 101, at z — 0. 
At redshift z — the simulation box contains 117 clusters 



with mass M(< i? v 
sample. 



r ) > 10 14 /i _1 M Q . This cluster set is our 



3.2 Simulated maps 

Around each cluster we extract a spherical region extending 
out to 6 7?vir. We then create maps of the relevant quantities 
along three orthogonal directions to investigate projection 
effects. Each map is a regular N p x N p grid. In the Tree+SPH 
code, each gas particle has a smoothing length hi and the 
thermodynamical quantities it carries are distributed within 
the sphere of radius hi according to the compact kernel 




< X SC i 

5 < a; < 1 
x > 1 , 



where x — r/hi and r is the distance from the particle posi- 
tion. We therefore distribute the quantity of each particle on 
the grid points within the circle of radius hi centered on the 
particle. Specifically, we compute a generic quantity on 
the grid point {j, k} as q jk dl = J q(r)dld 2 p = J2 Qi(mi/ Pi)Wi 
where d p is the pixel area, the sum runs over all the particles, 
and Wi oc J W(x)dl is the weight proportional to the frac- 
tion of the particle proper volume rrii/pi which contributes 
to the grid point {j, fc}. For each particle, the weights Wk are 
normalized to satisfy the relation ^ w k = 1 where the sum is 
now over the grid points within the particle circle. When hi 
is so small that the circle contains no grid point, the particle 
quantity is fully assigned to the closest grid point. Results 
in this paper are shown for a number of pixels Np = 128, 
corresponding to a comoving length resolution ~ 43/i _1 kpc, 
on average. Figures and [5] show an example of the thermal 
and the kinetic SZ maps of a cluster in our simulation. 



4 SCALING RELATIONS 

Galaxy clusters are self-similar in their dark matter com- 
ponent, but not necessarily in their baryonic component. 
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Non-gravitational processes, namely gas heating and cool- 
ing, vary the entropy of the ICM and can substantially affec t 
the distribution of the gas l|Borgani et alJl2004t IVoidl2004l . 
Measuring scaling relations among gas properties is an im- 
portant step to understand the role of the different processes 
governing the ICM physics. The list of the relevant processes 
that we included in our simulation is currently one of the 
most complete. Therefore, we can attempt a comparison of 
our simulated gas properties with the observations of the SZ 
effect. 

Here, we will investigate the relations Lx—yo and T — yo 
between the X-ray luminosity Lx or the ICM temperature 
T and the central Comptonization parameter yo; we will also 
compute the relation T — AS, where AS is the SZ surface 
brightness integrated over a given solid angle. 

The self-similar model of cluster formation predicts 
power-law scaling rela t ions between thes e quantities (e.gj . 
Evrard fe Henrvl Il99lt iKaiser 199 J Il99lt lEke et al.1 Il99a 
da Silva et al The cluster virial mass scales as M V i T oc 



p c (z)A c (z) R 3 , where R is the cluster radius in physical units 
and p c (z) is the critical density of the universe; pc(z) scales 
with redshift z as p c {z) oc E 2 (z) — fio(l + z) 3 + Q.a, in a flat 
universe. If we assume that clusters have a constant gas mass 
fraction and A c (z) is constant with redshift, 1 the cluster size 
R scales with z and the ICM mass M as Rx M 1/3 E~ 2/3 (z). 
Therefore the ICM mass scales with its temperature T as 
ME(z) oc T 3/2 ; for the CMB flux variation due to the ther- 
mal SZ effect we have AS oc J y(0)dtt oc d A 2 f Tn e d 3 r oc 
d A 2 T 5 ^ 2 E~ 1 (z), where d\ is the angular diameter distance 
to the cluster, and d 3 r = d\&Q,Al. We can also write AS 
in a different way to get its explicit dependence on yo: 
AS oc y d A 2 f dA oc y d- A 2 M 2/3 E- i/3 (z) oc y d A 2 TE- 2 (z), 
where dA = d 2 A dQ. is the infinitesimal projected area covered 
by the cluster. We find yo oc T 3/2 E(z). Finally, the X-ray 
luminosity Lx oc J n e riH A(T)dV oc T 2 E(z), if the cooling 



function A(T) oc T 1/2 , 
summary, we expect 



and we find yo oc L 



3/4 



E^ 4 (2 



r3/4 t-, 

yo oc Ly E 



1/4 



yo oc T 3/2 E(z) , 
ASd 2 A E(z) cc T 



i/2 



111 



(23) 
(24) 
(25) 



if the self-similar model holds. 

Figure|3]shows the correlation between the central peak 
yo of the Comptonization parameter map and the X-ray 
bolometric luminosity Lx- Since the simulation only in- 
cludes an approximate treatment of metal production and 
the cooling function used in the code is c omputed for zero 
metallicity, we follow llBorgani et al] i2004) and consider an 
optically thin gas of primordial cosmic abundance X = 0.76 
and Y = 0.24 to compute the ion number density nm+niiei- 
The X-ray luminosity of a simulated cluster is thus 



E 



n ei (riHi + nHei)A(Ti)dVi 



(26) 



In a flat universe, when f2o = 0.1, A c (z) increases by a factor 
~ 2.5 between, e.g., 2 = and 2 = 5, and this factor decreases 
with increasing Qo- On the other hand, E 2 (z) increases by a factor 
~ 22 between the same redshifts, and this factor increases with 
increasing Qo- 
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Figure 3. Central peak of the Comptonization parameter map 
vs. the X-ray bolometric luminosity. The dots are the simulated 
cluste rs, the open and solid c ircles wit h error bars th e cluster sam- 
ple of lMcCarthv et al] i2003rjl and ICooravl Il999fl . respectively. 
The solid line is the best fit to the simulated clusters and the 
two dotted lines show the ±3<r range from the best fit. The two 
dashed lines are the best fits to the ob served samples. The steep er 
dashed line is the fit to the sample of lMcCarthv et al] J2003b). 



with the appropriate cooling function A(T) 
iSutherland fc Dopital Il993tl : dVi = mi/ pi is the vol- 
ume occupied by the gas particle of mass rrii, density p; 
and electron number density n ei , Ti is its temperature 
and the sum runs over the gas particles, within R v i r , with 
Ti > 3 x 10 4 K and Pi smal l er than 500 times the mea n 
baryon density l|Croft et al.l l200ll iBorgani et ahl l2004l l. 
The sum therefore does not include dense star-forming 
particles, which are regulated by the multiphase model and 
may produce, when included, a spurious, highly uncertain 
contribution to the X-ray emission. 

The best fit to the simulated cluster sample (the dots 
in Figure |2J is 



log 



yo 



= (-0.57 ±0.01)+ 
+ (0.79 ± 0.02) log 



10 44 /i" 2 erg s" 1 



,(27) 



in agreement with the self-similar model. We do not attempt 
to simulate observational error bars on the cluster data ex- 
tracted from the simulation, because such a detailed analy- 
sis is beyond the scope of this paper. Therefore, in the ab- 
sence of errors on the simulated data, the uncertainties on 
the fitting param eters we compute are only indicative (e.g. 
iPress et al.lll992l p. 661). An alternative estimate of the un- 
certainty on this scaling relation is provided by the relative 
variance on a logarithmic scale a 2 — ([log(yo) — log^o')] 2 }. 
We find a = 0.10, comparable to the typical uncertainties 
of observed yo's. The dotted lines in Figure |S| show the ±3er 
amplitude. 

An analogous fit to 17 observed clusters (open 
circles in Figure 13 draw n from the sample com- 
piled by iMcCarthv et al] <2003bft. mostly based on 
the sample of Ir^ese^t all 120021 . yields the re- 
lation: log(yo/10 -4 ) = (-0.37 ± 0.05) + (0.65 ± 
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Figure 4. Central peak of the Comptonization parameter map vs. 
the X-ray emission weighted te mperature. Th e steep er dashed line 
is the best fit to the sample of lBenson et alJ l2004t) (solid circles 
with error bars). Other symbols and lines are as in FigureEl The 
temperatures of the real clusters are corrected for the presence of 
a cooling flow. 

0.04)log[£ 1/3 (z)Lx/10 44 /i~ 2 erg s" 1 ], in reasonable agree- 
ment with the self-similar model and our simulated sample. 
In computing E(z), we assume Qo — 0.3 and J1a = 0.7. 

In passing, we note that the apparent bimodal distribu- 
tion between simulated and observed data is a consequence 
of the observational bias of pointing at the very massive clus- 
ters: the volume of our simulation box is not large enough 
to include that many massive clusters. 

Despite their small statistical errors, the fit parame- 
ters on the real sa mple shou l d be treated with caution. In 
fact, the sample of ICooravl jl999l) (solid circles) yields a 
shallower relation log(t/ /10" 4 ) = (-0.08 ± 0.05) + (0.47 ± 
0.07)log[£ 1/3 (z)Lx/10 44 h- 2 erg s" 1 ]. 2 The discrepancy be- 
tween the two observed samples might be statistically mean- 
ingless, however, because the probability q of getting a larger 
X 2 with a sample drawn from a cluster population which 
obeys either fit to the observed samples is < 4 x 10 -8 (Table 
0. The small uncertainties claimed for the observed values 
explain the large x 2 ' s ; which on turn originate the small sig- 
nificance levels q. The discovery of possible systematic errors 
can of course bring the observed samples in better agreement 
with each othe r and with a single scaling relation. Indeed the 
ICooravl lll999l) sample has 10 out of 14 clusters in common 
with lMcCarthv et alJ J2003bl) sample, but this latter sample 
have more recent and different measures of yo and, in some 
cases, of Lx'- these differences might reflect both random 
and systematic errors. 

An intrinsic scatter in the yo — Lx scaling relation is 
in any case expected, due to the simplifying assumptions 
used to derive yo, namely the isothermality of the ICM, the 
f3— model and the cluster sphericity. It is remarkable that 
the scatters shown by our simulated clusters and by the 
observed samples are comparable. 



Figure 5. SZ flux decrement at 145 GHz, integrated within r^soOi 
vs. the X-ray emission weighted temperature. Symbols and the 
dashed line are as in Figure [I] The solid lines are the fits to 
the simulated sample using the least absolute deviation method. 
Dot-dashed lines show the range ±3D where D is the mean abso- 
lute deviation, from the fitting relation, of each simulated cluster 
in the AS-direction. The shallower relation, which corresponds to 
the shaded area, assumes Tx as the independent variable, whereas 
the steeper relation assumes E{z)ASd' A as the independent vari- 
able. The ±3-D range of the shallower relation has been shaded 
for clarity. The temperatures of the real clusters are not corrected 
for the presence of a cooling flow. 



Table 1. Fit parameters to the scaling relations log(/2V) 


= a + 


61og(/iX). 












X-Y 




a 


b 




Q 


Lx - yo 


(1) 


-0.08 ± 0.05 


0.47 ±0.07 


4.09 


2.5e-8 


fl = E 


(2) 


-0.37 ±0.05 


0.65 ±0.04 


4.31 


3.9e-8 


/2 = 1 


(4) 


-0.57 ±0.01 


0.79 ±0.02 






Tx - yo 


(1) 


-1.60 ±0.31 


1.87 ±0.31 


1.78 


2.5e-2 


fi = £ 2 


(2) 


-1.88 ± 0.39 


2.24 ±0.39 


1.06 


0.39 


/2 = 1 


(3) 


-2.31 ± 0.54 


2.79 ±0.51 


3.03 


1.7e-4 




(3a) 


-2.60 ±0.59 


3.27 ±0.60 


6.94 


1.3e-13 




(4) 


-1.31 ± 0.01 


1.55 ±0.03 








(3) 


2.73 ±0.38 


2.26 ±0.38 


0.70 


0.76 


h = 1 


(3a) 


2.07 ±0.42 


3.02 ± 0.44 


2.30 


4.9e-3 


h = d 2 A E 


(4) 


1.42 ±0.04 


2.41 ±0.11 








(4a) 


0.88 ±0.12 


3.86 ±0.19 







(1) ICooravl |l999l) : (2) iMcCarthv et alJ i2003rl : (3)/(3a) 
iBenson et alJ i2004l) with temperatures corrected/not corrected 
for the presence of a cooling flow; (4) simulated clusters; in the 
Tx — AS relation the fits to the simulated sample are computed 
with the method of the least mean absolute deviation; (4a) is 
the fit to the simulated sample assuming E{z)ASd? A , rather than 
Tx, as the independent variable. Units: Lx- 10 44 h~ 2 erg s — 1 ; Tx'- 
keV; yo'- 10 -4 ; ASd 2 A : Jy Mpc 2 . v is the number of degrees of 
freedom; q is the significance level of X* /v. 



2 Wherever appropriate, in compu ting the best fit we consider 
the uncertainties on both axes (e.g. iPress et al1ll992l p. 666). 
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Discrepancies among observed samples can also be seen 
in the relation between yo and the X-ray temperature Tx 
(Figure For the simulated clusters we plot the emission- 
weighted temperature averaged over the gas particles within 

Rvii 



T x 



X>ei("-Hi + nu a i)A(Ti)T i dV l 
J2 n ei( n m + nHa)A(Ti)dVi 
The simulated clusters yield the relation 



log( 



yo 

1Q- 1 



(-1.31±0.01) + (1.55±0.03) log 



keV 



(28) 



-(29) 



as expected in the self-similar model. For this relation, the 
standard deviation on a logarithmic scale is a = ([log(j/o) — 
l°g(yo t )} 2 ) 1/2 = 0.10 (the dotted lines in Figure |I] show the 
±3cr range). 

The sample of ICooravl il999T) yields a relation close to 
the self-similar prediction: log(j/o/10~ ) = ( — 1.60 ± 0.31) + 
(1.87 ± 0.31)log [ g 2/3 (z) Tx/keV1, whereas the sample of 
McCarthy et all J2003bti and the more recent sample of 
Benson et al.l l|2004 [). which ha s 12 out of 15 clusters in com- 
mon with the lMcCarthv et al. (2003b) sample, yield consid- 
erably steeper relations (Table [J. 

Some observed clusters have a central cooling flow. The 
temperatures of real clusters shown in Figure[I]are corrected 
for this effect. Neglecting this correction generally decreases 
the estimate of Tx and can substantially change the fit pa- 
rameters (Table 0. From the statistical point of view, how- 
ever, we can again draw ve ry little conclusions from most of 
these Tx — yo fits: only the lMcCarthv et alJ l)2003rJ) sample 
yields a large significance level q = 0.39. 

As we mentioned earlier, the discrepancies that we find 
among the different observed samples may be due to system- 
atic er rors affecting the data analysis. Indeed. iBenson et alJ 
i2004h have recently shown how the central yo's calculated 
from the measurements made with the SuZIE II receiver can 
be up to 60 per cent higher than the values de rived from the 
BIMA and OVRO interferometers. Therefore. IBenson et alJ 
(2004) prefer to use the SZ flux decrement, rather than the 
central yo, as an indicator of the magnitude of the SZ effect. 

Figure|^]shows the relation Tx— AS, where AS is the SZ 
flux within the circle of radius r2soo centered on the cluster: 
the average mass density within the sphere of radius r25oo is 
2500 t imes the critical density. According to IBenson et alJ 
(2004), this radius is a suitable choice for the instrumen- 
tal properties of Su ZIE II; moreover, r25oo is used in X-ray 
analyses of clusters ijEvrard. Metzler fc Navarrdllifflfih . The 
simulated sample shows a rather sparse distribution in the 
Tx — AS plane. Therefore we use the least absolute mean 
deviation method, which, in this case, is more appropriate 
than the least square fit method, to compute the fitting re- 
lation. We find 



log 



ASd 2 A 



Jy Mpc 2 



= (1.42±0.04) + (2.41±0.11) log , (30) 



as expected in the self-similar model. The observations agree 
with the self-similar slope only when we use the cooling flow 
corrected temperatures (sample 3 in Table Q. In this case, 
however, the normalization of the observed sample is sub- 
stantially larger than the model. Using the temperatures not 
corrected for the presence of the cooling flow, as shown in 
Figure (sample 3a) , yields a lower normalization but a 



substantially different slope. We do not have, however, error 
bars on our simulated data. Therefore, the fit is sensitive to 
the choice of the independent variable. Assuming E(z)ASd A 
as the independent variable in the simulated sample (sample 
4a) does indeed yield a substantially steeper slope; moreover, 
the sample variance is large enough to compensate the factor 
of two difference in the normalization. 

In conclusion, our simulated cluster sample seems to 
provide a reasonable description of the currently observed 
scaling relations between SZ and X-ray properties. To derive 
stronger constraints on our simulations of the ICM, however, 
we should perform more realistic mock observations of SZ 
clusters, although the disagreement, that we find among the 
observed samples of SZ clusters, casts some doubts on the 
robustness of the observed relations available to date. 



5 THE BASIC ASSUMPTIONS ON THE 
VELOCITY MEASUREMENTS 

The measure of the line-of-sight peculiar velocity of clusters 
with the SZ effects relies on two fundamental assumptions 
about the physical properties of the ICM: (1) bulk flows of 
gas clouds within the ICM are negligible and gas and dark 
matter have the same mean peculiar velocity, and (2) the 
ICM temperature estimated from the X-ray spectrum is a 
reliable estimate of the temperature weighted by the elec- 
tron number density. We now check on these assumptions 
on turn. 



5.1 Velocity bias and internal bulk flows 

We might expect that turbulence induced by mergers, 
gas cooling, galactic winds and stellar feedback can pro- 
vide substantial contribution to the bulk flow of gas 
clouds within the ICM (e.g. [i nogamov fc Sunv aevl 120031 : 
ISunvaev. Norman fc Brvanll2003l: ISchuecker et alJl2004|) ~ 
clusters with low peculiar velocity, these effects can also re- 
sult in a gas mean velocity quite different from the mean 
velocity of the dark matter halo. Our simulation, which in- 
cludes these physical processes, should provide a sample of 
realistic clusters for the investigation of these effects. 

In terms of velocity magnitudes, for each simulated clus- 
ter, the three components of the mean velocity of the gas 



particles within R vl[ 



never deviate substantially from 



the corresponding components of the mean velocity of the 
dark matter particles, ddm: we find a mean absolute devia- 
tion, a more robust estimator of the width of a distribution 
than the standard deviation, (|w gas — «dm|) = 18±15kms _1 ; 
the maximum value in our sample is |w gas — «DM|max = 88 
km s . 

However, in terms of relative velocities, the bias can be 
substantial. Setting a threshold at i>fj M = 100 km s _1 , we 
find a mean deviation {5} = (|(%as — i>dm)/vdm|) = 0.9±2.2 
for the clusters with ddm < ^dm, an order of magni- 
tude larger and with a substantially larger deviation than 
(5) = 0.09 ± 0.08 for the clusters with wdm > i>dm- We con- 
clude that we should not expect any relevant systematics 
from velocity bias when we consider clusters with line-of- 
sight peculiar velocities larger than ~ 100 km s _1 . 

We also need to quantify how much the bulk flow of 
gas clouds within the ICM can affect the peculiar velocity 
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Figure 6. Average (vij) of the velocities at the pixel sites vs. the 
peculiar velocity of the cluster ^dm- Only pixels within -R v i r of 
the cluster centre are included in the calculation of the means. 
The dotted line is the relation (vij) = ^dm- 
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Figure 8. Electron density weighted temperature T e vs. X-ray 
emission-weighted temperature Tx- Only gas particles within 
ri; m of the cluster centre are considered. Each panel shows the 
Tx — T e relation for a different value of ri; m , as indicated. The 
dotted line is the equation T e = Tx; the solid line is the best fit 
T e = a + bTx ■ The parameter a is in keV. 
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Figure 7. Distribution of the mean absolute deviations between 
the velocities at the pixel sites and the peculiar velocity of the 
cluster v = «dm (solid histogram) or the mean velocity of the 
velocity map v = (dashed histogram). Only pixels within 

Rvir of the cluster centre are included in the calculation of the 



estimate. We consider the velocity maps Vij = cfiij, where 
Pij is a generic element of the array corresponding to the 
map computed according to equation 1101 : therefore, the 
velocity Vij is a velocity weighted by the electron number 
density. We first compute the mean (vij) of the Vij values on 
the pixels within i? v i r of the cluster centre. This mean agrees 
with the true vdm (Figure [SJl: the mean absolute deviation 
is (\{vij) — «dm[) = 31 ± 32 km s _1 . 

For each map, we then compute the mean absolute de- 
viation {\vij — v\) of the velocity values at the pixel sites 
from v = vdm or v = (vij). Figure Q shows the distribution 
of the (\vij — v\)'s for all the maps. The average values are 
(\vij-v D M\) = 119T48 kms" 1 and = 115T45 



km s _1 . We do not expect the two distributions to be dif- 
ferent because (|(My) — vdm|) is small. The two samples, 
with wdm smaller or larger than the threshold 100 km s _1 , 
have a {\wj — v\) distribution similar to the complete sam- 
ple, both when v = udm and v — (vij). This last result 
indicates that the velocity bias found above for the slow 
clusters ({5) = 0.9 ± 2.2) is totally due to the bulk flow of 
gas clouds within the ICM, and not to different mean veloc- 
ities of gas and dark matter particles. This conclusion also 
explains why the scatter of (5) for the slow clusters (2.2) 
is much larger than the scatter for the fast clusters (0.08): 
bulk flows have roughly the same magnitude, independently 
of the cluster speed. 

According to FigureQ internal bulk flows produce a de- 
viation, from vdm, smaller than 200 km s _1 in 93 per cent 
of the cases. We conclude that the velocity bias and inter- 
nal bulk flows set a conservative lower limit of ~ 200 km 
s _1 to the unc ertainty on peculiar velocity m easures with 
the SZ effects jNaeai et alJl2003t lHoIderll2004l) . well below 
the uncertainties of th e currently most accu r ate observations 
j jLaRoaue et al j 120031 : iBenson et alJ l2003t iKitavama et alJ 
l2004F ""but similar to those expected in future measurements 
at microkelvin sen sitivity and arcminute angular resolution 
jKnox et al.ll2004) . 



5.2 Electronic and X-ray temperatures 

For all X-ray emitting particles, we can compute an X-ray 
emission weighted temperature (equation 1281 and an elec- 
tron density weighted temperature 



T e = 



n ei TidVi 



(31) 



In both equations 1281 and (131ft . the sums run only over the 
particles within a maximum radius rii m - 
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Figure 9. Temperature profiles of the simulated clusters. The 
upper left (right) panel shows the temperature weighted by the X- 
ray emission (electron number density). The bottom panel shows 
the profile of the ratio of the two temperatures. Dotted lines are 
the individual profiles. Dots with error bars show the mean pro- 
files with 1-a standard deviations. 



Note that T e equals a mass- weighted temperature Tm = 
Y^miTi/Y^ 171 ! when the gas is fully ionized, namely n e i oc 
pi, which is usually the case for temperatures larger than 
0.1 keV. 

Figure |H] shows T e vs. Tx for different values of 
rn m . The two temperatures are comparable only at the 
smallest ri j m , when T e is roughly 10 per cent larger 
than T x jMathiesen fc Evrardl boOll) . At larger nim's 
the disagreement increases. In fact, the ICM tempera- 
ture profile drops at large radii both in real clu sters 
(e.g. iDe Grandi fc Molendal20ol iKaastra et alJl2005t and 
in simulations (e.g. | Lm ^tai Tl200^ r TBor^anie^aiTl2004l: 
iRasia. Tormen fc Moscardinill2004) . 

Figure shows the profiles of our simulated clusters. 
Tx is largely insensitive to rum because it is proportional 
to the square of the gas density and is therefore dominated 
by the emission from the central region. Instead, T e steeply 
decreases with increasing ru m because it is only proportional 
to the first power of the electronic gas density. 

When rii m = i? v ir, T e can be substantially smaller than 
Tx'- the mean relation between these two measures of the 
ICM temperature is T e /keV = 0.17 + 0.69(T x /keV) (Figure 
[SJ. Therefore, in spatially poorly resolved clusters, using Tx 
rather than T e can substantially overestimate the peculiar 
velocity (equations 1151 and 1221 , as we will see in the next 
section. 

Note th at here we compute th e bolometric temperature. 
Figure 3 in iBoreani et al.l (|2004 ) shows that the tempera- 
tures estimated by weighting the emissivity in a finite en- 
ergy band are larger than those obtained by weighting with 
the bolometric emissivity. Therefore, for real clusters, the 
T e — Tx relation can be shallower and the overestimate of 
the peculiar velocity more severe. It remains to quantify this 
trend when a realis tic observational proc edure is applied to 
simulated clusters iMazzotta et alJl2004l) . 



6 PECULIAR VELOCITY MEASUREMENTS 

Measuring the peculiar velocity of clusters with the SZ ef- 
fect requires the separation of the thermal from the kinetic 
contribution to the total SZ flux. Spatially resolved and un- 
resolved clusters require two different procedures for this 
separation. 

In the former, a model for the ICM distribution is de- 
rived from the X-ray surface brightness to estimate the cen- 
tral values of the SZ contributions. In the latter, the beam- 
averaged SZ contributions are measured without any infor- 
mation on the ICM distribution. To date, only the former 
case has been applied to observations. 

We now consider these two approaches on turn to es- 
timate the systematics that they introduce on the peculiar 
velocity measurement. Here, we will only consider the phys- 
ical properties of the ICM and we will not model the instru- 
mental or data reduction problems which affect real obser- 
vations. 



6.1 Resolved clusters: fitting profile method 

We perform a rather simple simulation of the observation 
of individual clusters. Namely, we (1) construct an X-ray 
surface brightness map from the simulation, (2) derive the 
gas distribution shape parameters j3 c and 6 C from the X-ray 
surface brightness profile, (3) compute the SZ fluctuations in 
the CMB intensity maps AI(u) at four different frequencies 
v — 30, 143, 217, and 353 GHz, (4) derive the central values 
AJo(^) of the radial fits with the shape parameters /3 C and 8 C 
imposed by the X-ray surface brightness profile, (5) perform 
a fit to equation 112H by assuming that the four AIo(v)'s 
are independent, (6) estimate the peculiar velocity (equation 
I15H by using the X-ray temperature Tx (equation or the 
electronic temperature T e (equation 1311 . 

In equation 1151 . T e indicates the central ICM temper- 
ature weighted by the electron density. Observationally, one 
uses the X-ray temperature averaged within the resolution 
beam. To quantify the effect of the beam size, we compute 
the temperature with the gas particles within spheres of dif- 
ferent radii ta- 

Figure 1101 shows the distribution of the deviations 
(n os t — WDiiJ/fDM of the estimated line-of-sight peculiar ve- 
locity ii ost from the actual cluster velocity vdm- The solid 
and dashed histograms show these deviations when we use 
Tx or T e , respectively. We consider the temperatures aver- 
aged within four different radii. 

The temperature within the smallest radius 0.1i? v ir is 
the closest to the central temperature that should be used in 
equation 115> . At this radius, both the X-ray and the elec- 
tronic temperature overestimate the peculiar velocity. The 
overestimate with T e is more severe, because T e is 10 per 
cent larger than Tx (Figure [5J . The incorrect estimate of 
the peculiar velocity derives from the fact that the gas is 
neither spherically and smoothly distributed nor isothermal, 
as assumed when we impose the /3-model to the X-ray sur- 
face brightness. In fact, the agreement between estimated 
and real peculiar velocity does not improve when we use the 
/3-model parameters derived from the SZ surface brightness, 
rather than the X-ray surface brightness, despite the fact 
that the two parameter sets do not generally coincide (e.g. 
lYoshikawa. Itoh fc Suto1ll99StlLin et alJl2003l) . 
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Figure 10. Distributions of the deviations of the estimated line-of-sight velocities from the real velocities of spatially resolved clusters 
with VuM > 100 km s . The velocities are estimated with equations 1151 where T e = (Tx)r A (solid histograms), or T e = (T e ) rA (dashed 
histograms) and the /3-model parameters are derived from the fit to the X-ray surface brightness. Each panel shows the distributions for 
a different value of ta, as indicated. The other numbers in each panel are the median and the mean absolute deviation from the median 
of each distribution. 



The radial profile of the X-ray temperature Tx de- 
creases less rapidly than the radial profile of T e (Figure 
Therefore, the peculiar velocity estimate based on Tx 
is rather insensitive to the beam size, whereas the velocity 
estimated with T e rapidly decreases with ta- 

The distributions in Figure [T01 only include the clusters 
with «dm > 100 km s _1 . Including all the clusters leaves 
the medians unchanged, but substantially extends the tails 
of the distributions, because, in slow clusters, errors of order 
a few tens of km s _1 produce large relative deviations due 
to the bulk flows internal to the ICM. 

In terms of velocity magnitudes, the deviations are com- 
parable, on average, to the systematics due to the internal 
bulk flows. The upper panels in Figures [T?l show the scat- 
ter plot between the estimated and the real velocities for 
all the clusters when ta = 0.1i? v ir- The standard deviations 
are a — ((v cst — udm) 2 ) 1 ' 2 ~ 153 km s _1 and 173 km s _1 , 
when we use the X-ray or the electronic temperatures, re- 
spectively. 

If our simulated clusters mirror a realistic sample, we 
expect that, for spatially resolved clusters with peculiar ve- 
locity larger than ~ 100 km s _1 , the usual assumptions 
on the ICM physical properties should cause on average a 
~ 10 — 20 per cent systematic overestimate of the veloc- 



ity. For these clusters, the width of the distribution in the 
upper-left panel of Figure I1UI suggests ~ 50 per cent as a 
typical uncertainty on the velocity estimate. 

6.2 Unresolved clusters: aperture method 

When the radio/microwave observations do not spatially re- 
solve the cluster, we need to use the equations described in 
Sect. 12.21 Here, we compute the y(9) and b(0) maps, and 
compute AS(x) by integrating AI(x) over a beam of phys- 
ical size 2rA(@b) ~ dA0b- We then extract Y{6b) and B(0b) 
from the spectral fit. We emphasize that in this case there 
is no need of any assumption about the ICM spatial distri- 
bution. 

Figure 1111 shows the distributions of the deviations of 
the estimated line-of-sight velocities from the real velocities 
for four different aperture radii ta- Our maps have 128 x 128 
pixels and are 12i? v i r on a side. Therefore, we choose ta 
0.3-Rvir to avoid the shot-noise due to the small number of 
pixels in smaller apertures. 

When we use the X-ray temperature averaged over the 
appropriate aperture in equation I22II . the larger is the aper- 
ture, the larger is the line-of-sight peculiar velocity overes- 
timate (solid histograms). In fact, the ratio of the two inte- 



Velocities with the Sunyaev-Zeldovich effects: scaling relations and systematics 11 



- 1 — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — r 

— ^™ r A =d.3R; ir 



t — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — r; 

medL dev. r A =0.5R„w. ; 

. 0.24 0.23 I _ j 

0.12 0.19 " ' 

I i 
i i 



0.3 ^ — 



0.2 



0.1 
0.0 

0.3 

0.2 

0.1 
0.0 



med. dev. 
0.20 0.30 
. 0.14 0.29 



med. dev. 
_ 0.38 0.28 
. 0.07 0.16 , 



r A =1.0R vir 




med. dev. — 3.0R v j r 

0.52 0.44 




-1.0-0.5 0.0 0.5 1.0 -1.0-0.5 0.0 0.5 1.0 
(estimated vel. v est )/(actual vel. v DM )-l 



Figure 11. Distributions of the deviations of the estimated line-of-sight velocities from the real velocities of spatially unresolved clusters 
with vdm > 100 km s . The velocities are estimated with equations H22I where {T e )g b = (Tx)r A (solid histograms), or (T e )g b = {T e ) rA 
(dashed histograms). Each panel shows the distributions for a different value of r A , as indicated. The other numbers in each panel are 
the median and the mean absolute deviation from the median of each distribution. 



grated SZ components B(6b)/Y(0b), in units of the central 
ratio 60/2/0, increases by 45 per cent between ta = 0.3 and 
r.4 = 3-Rvir, while Tx remains basically constant. On the 
other hand, the rapid decrease of the electron temperature 
T e with ta compensates the increase of B(8b) /Y (6b), and 
the peculiar velocity estimate is basically unbiased (dashed 
histogram) . 

As for the case of spatially resolved clusters, Figure fTTl 
only shows the clusters with voai > 100 km s _1 : slower 
clusters produce longer tails in the distributions but do not 
vary the central values. 

The lower panels in Figure ED show the extreme case of 
unresolved clusters, with beam size ta = 3J? v i r . Even in this 
case the standard deviation is a — ((v cst — fDivi) 2 ) 1 ^ 2 < 160 
km s _1 and therefore comparable to the systematic errors 
due to the internal bulk flows. 

Figure fTTl clearly shows that, with only the X-ray tem- 
perature available, it is essential to have good spatial reso- 
lution, in order to have a measure of the ratio B(8b)/Y(9b) 
within a sufficiently small area where Tx and T e are com- 
parable. Otherwise the large value of Tx can severely over- 
estimate the peculiar velocity. 



7 CONCLUSIONS 

In principle, the evolution of the velocity field on large scales 
can be measured by detecting the kinetic SZ effect of galaxy 
clusters. There exist substantial technical difficulties in this 
measurement due to the weakness of the signal and due to 
the fact that the kinetic effect does not produce a spectral 
distortion of the CM B. Nevertheless, a few detections have 
been claimed to dat e llHolzapfel et alJ|l997l;|L aRoaue et alJ 
l2003t iBenson et aljl2003t iKitavama et alJl2004l) . and more 
sensitive instruments will make this measurement less prob- 
lematic in the near future. 

The estimate of the peculiar velocity of a cluster from 
the amplitude of the kinetic SZ effect assumes an isothermal 
ICM in hydrostatic equilibrium i n the cluster potentia l well. 
Recent X-ray observations (e.g. iForman et al.ll2003tl show 
that the ICM is actually far from this ideal model. 

We extracted a sample of 117 fairly realistic galaxy clus- 
ters, with mass larger than 1O 14 /i _1 M0, from an iV-body 
hydrodynamical simulation of a large volume of a ACDM 
cosmological model to estimate the systematic errors, due 
to the gas bulk flows and ICM temperature estimate, which 
can affect the measures of the peculiar velocity. 

Our sim ulation reproduces the observed proper ties of X- 
ray clusters iBorgani et alJl2004l : lEttori et al.ll2004|) and the 
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Figure 12. Relation between the estimated and the actual ve- 
locity. Upper (lower) panels show spatially resolved (unresolved) 
clusters. In left (right) panels the ICM temperature is weighted 
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where a 2 = ((»dm — i>est) 2 ). The numbers in each panels are the 
standard deviations a in km s _1 . 



distribution of the intracluster light jMurante et al.l I2004T) 
reasonably well. Here, we show that our model also repro- 
duces the observed scaling relation between X-ray and SZ 
properties, namely the X-ray luminosity or the emission- 
weighted X-ray temperature vs. the central value j/o of the 
Comptonization parameter map. The results of the com- 
parison, between observed and simulated clusters, that we 
propose here, are not yet fully satisfactory, mostly because 
of the lack of a uniform sample of clusters observed both in 
the microwave and X-ray bands. 

We then show that if the ICM is not in hydrostatic equi- 
librium in the cluster potential well, bulk flows of gas within 
the ICM itself can introduce a systematic error in the esti- 
mate of the peculiar velocity. However, this effect is unlikely 
to be larger than ~ 200 km s _1 . Our result, based on a 
large sample of simulated clusters, provides a robust st atis- 
tical confirma tion of the analyses of iNaeai et alJ i2005tl and 
iHolderl (I2004T) who considered the iV-body hydrodynamical 
simulations of only one and three clusters, respectively. 

We finally investigate the relevance of the estimate of 
T e , the ICM temperature weighted by the electron number 
density, in the peculiar velocity measurement. A possibile 
systematic error originates from the fact that the observed 
temperature Tx , which is derived from the X-ray spectrum, 
might be a biased estimate of T e . In the simulations, Tx 
is usually identified with the temperature weighted by the 
X-ray emission. When averaged within large radii of the clus- 
ter centre, Tx can overestimate T e by ~ 20 — 50 per cent. If 
the angular resolution of the radio/microwave observation 
is too poor to provide a spatially resolved profile of the SZ 
surface brightness and it only provides an integrated flux, 
the temperature overestimate propagates into a peculiar ve- 
locity overestimate of the same amount. 

On the other hand, in the case of clusters which are 
spatially resolved in the microwave bands, we can use the 



temperature of the ICM in the central region of the clus- 
ter, where Tx is comparable to T e . For observations which 
resolve a tenth of the virial radius, the velocity is still over- 
estimated by 10 — 20 per cent, because of the assumptions 
on the spatial and thermal distribution of the ICM, usually 
parametrized with a /3-model, used to separate the thermal 
and kinetic SZ components. 

Our results rely on the usual assumption that the 
emission-weighted temperature extracted from the simu- 
lations is comparable to the ICM temperature derived 
from the observed X-ray spectrum. The emission-weighted 
temperature can however over estimate the spectral tem- 
perature by 20 - 30 per cent llMathiesen fc Evrardl l200ll : 
iMazzotta et alJl2004fc iRasia et alJl2004al) . This uncertainty 
makes the identification of the observed Tx with T e more 
problematic, although it seems to alleviate the systematics 
we find here. 

By combining all these results we can set a conserva- 
tive lower limit of ~ 200 km s _1 (see Figures |7| and 1121 
to the systematic errors that will affect upcoming measure- 
ments of the cluster peculiar velocity with the SZ effect. 
This limit is comparable to the accuracy expected when 
the contamination from dusty galaxies and primary CMB 
anisotropies is considered in radio/microwave observations 
at microkelvin sen sitivity and arcminute angular resolution 
jKnox et al.ll2004) . Therefore, the internal ICM bulk flows 
and its temperature estimate may limit any substantial im- 
provements that could come from more sensitive microwave 
measurements and more accurate component separations. 

We consider the systematic error we find here a lower 
limit because we identify clusters in the three-dimensional 
space of the simulation box and analyze their X-ray and SZ 
properties in two-dimensional projected maps without in- 
cluding instrumental noise and back-/fore-ground sources. 
Real-world complicati ons can add further systematic errors 
iAghanim et al.l|2004l) . More realistic mock observations of 
SZ clusters are required to design a method to correct for 
these biases. To accomplish this task one can construct mi- 
cro wave sky maps from iV-body hydrodynamical simulations 
fe.g.lSch afer et all2004l) . identify clusters in these maps (e.g. 
iMelin et alJl2004T) and analyze the simulated data according 
to standard observational procedures. This approach will 
help to quantify how, in a realistic situation, the different 
systematic errors propagate into the statistical analysis of 
cluster velocity fields. We plan to investigate this issue in 
future work. 



ACKNOWLEDGMENTS 

We thank Gil Holder and Ian McCarthy for pointing out 
a few inaccuracies, regarding the observational data, that 
were present in an early version of the paper, and the ref- 
eree, Douglas Scott, for relevant suggestions that improved 
the presentation of our results. The simulation ran on an 
IBM-SP4 machine at the "Consorzio Interuniversitario del 
Nord-Est per il CAlcolo elettronico" (CINECA, Bologna), 
with an INAF-CINECA CPU time grant. This work has 
been partially supported by the NATO Collaborative Link- 
age Grant PST.CLG. 976902, the INFN Grant PD-51, the 
MIUR Grant 2001, prot. 2001028932, "Clusters and groups 
of galaxies: the interplay of dark and baryonic matter", and 



Velocities with the Sunyaev-Zeldovich effects: scaling relations and systematics 13 



by ASI. KD acknowledges support by a Marie Curie Fellow- 
ship of the European Community program Human Potential 
under contract number MCFI-2001-01227. 



REFERENCES 

Aghanim N., Gorski K. M., Puget J.-L., 2001, A&A, 374, 
1 

Aghanim N., Hansen S. H., Lagache G., 2004, A&A, sub- 
mitted, astro-ph/0402571 

Atrio-Barandela F., Kashlinsky A., Miicket J. P., 2004, 
ApJ, 601, Llll 

Benson A. J., Reichardt C, Kamionkowski M., 2002, MN- 
RAS, 331, 71 

Benson B. A. et al., 2003, ApJ, 592, 674 

Benson B. A., Ade P. A. R., Bock J. J., Ganga K. M., 
Henson C. N., Thompson K. L., Church S. E., 2004, ApJ, 
submitted, astro-ph/0404391 

Borgani S. et al., 2004, MNRAS, 348, 1078 

Carlstrom J. E., Holder G. P., Reese E. D., 2002, ARA& 
A, 40, 643 

Cavaliere A., Menci N., 2001, MNRAS, 327, 488 

Cohn J. D., Kadota K., 2004, astro-ph/0409657 

Cooray A. R., 1999, MNRAS, 307, 841 

Courteau S., Dekel A., 2001, in Astrophysical Ages and 

Times Scales, ASP Conference Series Vol. 245. edited by 

T. von Hippel, C. Simpson and N. Manset, San Francisco: 

Astronomical Society of the Pacific, p. 584 
Croft R. A. C, Di Matteo T., Dave R., Hernquist L., Katz 

N., Fardal M. A., Weinberg D. H., 2001, ApJ, 557, 67 
da Silva A. C, Barbosa D., Liddle A. R., Thomas P. A., 

2001, MNRAS, 326, 155 
da Silva A. C, Kay S. T., Liddle A. R., Thomas P. A., 

2004, MNRAS, 348, 1401 
De Grandi S., Molendi S., 2002, ApJ, 567, 1 
Diaferio A., Sunyaev R. A., Nusser A., 2000, ApJ, 533, L71 
Diaferio A., Nusser A., Yoshida N., Sunyaev R. A., 2003, 

MNRAS, 338, 433 
Diego J. M., Mohr J., Silk J., Bryan G., 2003, MNRAS, 

341, 599 

Dos Santos S., Dore O., 2002, A&A, 383, 450 
Eke V. R, Cole S, Frenk C. S, 1996, MNRAS, 282, 263 
Eke V. R., Navarro J. F., Frenk C. S., 1998, ApJ, 503, 569 
Ettori S., et al., 2004, MNRAS, 354, 111 
Evrard A. E., Henry J. P., 1991, ApJ, 383, 95 
Evrard A. E., Metzler C. A., Navarro J. F., 1996, ApJ, 469, 
494 

Feldman H. A., et al. 2003, ApJ, 596, L131 

Forman W., et al., 2003, in High Energy Processes 

and Phenomena in Astrophysics, IAU Symposium 214, 

astro-ph/0301476 
Forni O., Aghanim N., 2004, A&A, 420, 49 
Haehnelt M. G., Tegmark M., 1996, MNRAS, 279, 545 
Hamilton A. J. S., 1998, in The Evolving Universe. Selected 

Topics on Large- Scale Structure and on the Properties of 

Galaxies, Dordrecht: Kluwer Academic Publishers, p. 185 
Hansen S. H., 2004, MNRAS, 351, L5 
Hansen S. H., Pastor S., Semikoz D. V., 2002, ApJ, 573, 

L69 

Holder G. P., 2004, ApJ, 602, 18 

Holder G. P., Haiman Z., Mohr J. J., 2001, ApJ, 560, Llll 



Holzapfel W. L., Ade P. A. R., Church S. E., Mauskopf P. 
D., Rephaeli Y., Wilbanks T. M., Lange A. E., 1997, ApJ, 
481, 35 

Huterer D., Kim A., Krauss L. M., Broderick T., 

astro-ph/0402002 
Inogamov N. A., Sunyaev R. A., 2003, Astron. Lett., 29, 

791 

Itoh N., Kohyama Y., Nozawa S., 1998, ApJ, 502, 7 
Kaastra J. S., et al., 2003, A&A, 413, 415 
Kaiser N., 1991, ApJ, 383, 104 

Kashlinsky A., Atrio-Barandela F., 2000, ApJ, 536, L67 
Kitayama T., et al, 2004, PASJ, 56, 17 
Knox L., Holder G. P., Church S. E., 2004, ApJ, 612, 96 
LaRoque S. J., Carlstrom J. E., Reese E. D., Holder G. P., 

Holzapfel W. L., Joy M., Grego L., 2003, ApJ, submitted, 

astro-ph/0204134 
Lin W. C, Nor man M. L., Bryan G. L., 2003, 

astro-ph/0303355 
Mather J. C, Fixsen D. J., Shafer R. A., Mosier C, Wilkin- 
son D. T., 1999, ApJ, 512, 511 
Mathiesen B. F., Evrard A. E., 2001, ApJ, 546, 100 
Matsubara T., Szalay A. S., 2002, ApJ, 574, 1 
Mazzotta P., Rasia E., Moscardini L., Tormen G., 2004, 

MNRAS, 354, 10 
McCarthy I. C, Babul A., Holder G. P., Balogh M. L., 

2003a, ApJ, 591, 515 
McCarthy I. C, Holder G. P., Babul A., Balogh M. L., 

2003b, ApJ, 591, 526 
Mei S., Bartlett J. G., 2003, A&A, 410, 767 
Mei S., Bartlett J. G., 2004, A&A, 425, 1 
Melin J.-B., Bart lett J. C, Delabrouille J., 2004, 

astro-ph/0409564 
Moscardini L., Bartelmann M., Matarrese S., Andreani P., 

2002, MNRAS, 335, 984 
Murante G., et al., 2004, ApJ, 607, L83 
Nagai D., Kravtsov A. V., Kosowsky A., 2003, ApJ, 587, 

524 

Pointecouteau E., Giard M., Barret D., 1998, A&A, 336, 
44 

Press W. H., Teukolsky S. A., Vetterling W. T., Flannery 
B. P., 1992, Numerical Recipes in C - The Art of Scientific 
Computing, second edition, Cambridge University Press 

Rasia E., Mazzotta P., Borgani S., Moscardini L., 
Dolag K., Tormen G., Diaferio A., Murante G., 2004a, 
astro-ph/0409650 

Rasia E., Tormen G., Moscardini L,, 2004b. MNRAS. 351, 
237 

Reese E. D., Carlstrom J. E., Joy M., Mohr J. J., Grego 

L., Holzapfel W. L., 2002, ApJ, 581, 53 
Rephaeli Y., 1995, ARA&A, 33, 541 

Schafer B. M., Pfrommer C, Bartelmann M., Springel V., 

Hernquist L., 2004, astro-ph/0407089 
Schafer B. M. , Pfrommer C, Zaroubi S., 2003, MNRAS, 

submitted, astro-ph/0310613 
Schuecker P., Finoguenov A., Miniati F., Bohringer H., 

Briel U. G., 2004, A&A, submitted, astro-ph/0404132 
Schulz A. E., White M., 2003, ApJ, 586, 723 
Sheth R.K., Diaferio A., 2001, MNRAS, 322, 901 
Springel V., Hernquist L., 2002, MNRAS, 333, 649 
Springel V., Hernquist L., 2003, MNRAS, 339, 312 
Springel V., Yoshida N., White S. D. M., 2001, NewA, 6, 

79 



14 Diaferio et al. 



Sunyaev R. A., Norman M. L., Bryan G. L., 2003, Astron. 
Lett., 29, 783 

Sunyaev R. A., Zeldovich Ya. B., 1972, Comments Astro- 

phys. Space Phys., 4, 173 
Sutherland R. S., Dopita M. A., 1993, ApJS, 88, 253 
Verde L., Haiman Z., Spergel D. N., 2002, ApJ, 581, 5 
Voit M., 2004, Rev. Mod. Phys., in press 

(astro-ph/0410173) 
Weller J., Battye R. A., Kneissl R., 2002, PhRvL, 88, 

231301 

White M., Hernquist L., Springel V., 2002, ApJ, 579, 16 
Yoshikawa K., Itoh M., Suto Y., 1998, PASJ, 50, 203 



This figure "fOl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405365v2 



This figure "f02.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0405365v2 



